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DESCRIPTIVE EXAMPLES OF GRAVITY ASSISTED TRAJECTORIES 

by 

James B. Eades, Jr.* 

Special Projects Branch 


ABSTRACT 

Two approximate methods for calculating the planetary Swing-By 
maneuver are presented. These schemes may be used to determine, 
analytically, the path of motion, time of motion, energy and velocity 
change as experienced by a satellite prior, during and often after the 
encounter. 

Also, an examination for the extremals of the energy and velocity 
changes is discussed as a part of the overall Swing-By problem. 


NAS-NASA Resident Research Associate, on leave from the Virginia Polytechnic 
Institute, Blacksburg, Va. 
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DESCRIPTIVE EXAMPLES OF GRAVITY ASSISTED TRAJECTORIES 


I. INTRODUCTION 

Gravity assisted trajectories, or ”swing-by” maneuvers, are produced by 
the passage of a satellite (man-made or otherwise) close to a massive celestial 
body. As a consequence of this the satellite’s flight path is perturbed, significantly, 
due to the mass attraction of the larger body. For situations such as these the 
flight path develops something akin to a ”kink;” the vehicle will have its energy 
of motion altered, and its velocity vector will undergo a change in direction. 

These several happenings are all due to the gravitational attraction of one moving 
massive body on a second smaller moving mass, (say) the satellite. 

As an example of the problem type which will be described herein, consider 
Figure 1.1 where an earth launched satellite is depicted to pass close to one of 
the massive outer planets. Because of the close proximity of the spacecraft to 
the attracting planet, the vehicle’s path is altered. 

It is not difficult to visualize that one might be able to ’’control” the effects 
of such a swing-by operation; especially if it would be possible to ’’control” how 
close the satellite approached to the planet, and if some ’’control” could be main- 
tained over the speed of the vehicle. As a matter of interest previous theoretical 
studies, related to space travel, have indicated that significant fuel savings and/or 
increased payloads - as well as shortened flight times - could be achieved by a 
proper utilization of this technique. 
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v d = HELIOCENTRIC VELOCITY OF 
SWING -BY PLANET (fT^). 

V_? HELIOCENTRIC VELOCITY OF 
S/C BEFORE SWING- BY. 


V+= S/C HELIOCENTRIC VELOCITY 
AFTER SWING- BY 

av l !s/c VELOCITY INCREMENT 
FOR LAUNCH. 

AVd= S/C VELOCITY INCREMENT 
DUE TO SWING- BY 



Figure 1.1-Sketch Depicting the Coplanar Circular Orbits of Earth and the Sw*ng-By 
Planet, with effects of Swing-By (AVj and Path Change) Illustrated. 


II. BASIC ASSUMPTIONS 


In order to simplify the mathematics for this investigation, and to aid in pro- 
viding analytic solutions, several assumptions will be made. In general, these 
are as follows: 

1. All motions will be treated as two-body situations; thus, when the satellite 
is in the vicinity of a disturbing planet it will be assumed that the planet provides 
the principal influence on the trajectory. Consequently, all other bodies are 
neglected and the problem is referred to planetocentric space. However, for 
motion in interplanetary space it will be assumed that the primary influence is the 
Sun and the study is (then) referred to heliocentric space. 

These assumptions naturally lead to the conceptual idea of a ’’Sphere of 
Influence,” which will be employed when the gravitational influence of an attracting 
body, on the flight path, is that of the disturbing mass. 

2. During that time period when the vehicle is in passage through the sphere 
of influence, the attracting planet is presumed to move only negligibly along its 
own orbit; hence, the ”swing-by” is presumed to occur at a fixed heliocentric 
location for the attracting body. 

3. The trajectory flown by the satellite through the Sphere of Influence will 
be a hyperbola. This will be treated as a constant energy path; thus the vehicle 
will be influenced solely by the attracting planet. For this condition, the trajectory 
will be described by its constant specific energy and a fixed specific angular 
momentum (relative to the attracting planet). 

NOTE: The interplanetary trajectories, under the same restrictions as noted 
above, are two-body flight paths, where the satellite moves about under the singular 
influence of the Sun. In this regard the vehicle’s trajectory, outside of the attract- 
ing planet's sphere of influence, is made up of heliocentric arcs characterized 
by constant energy and angular momentum (so long as the only force considered 
is that of the Sun's gravitational attraction.) 

It should be mentioned that the solution technique employed here is by no 
means exact; it is an approximation utilized to alleviate some of the computational 
difficulties which necessarily arise when the true multi-body aspects of this 
problem are introduced. The case studies included here are variations of the 
patched conic mode of solution; the two modes of approach which constitute this 
investigation are included to indicate the degrees of complexity which can be intro- 
duced into this simpler analytic approach. Needless to say, when extreme accuracy 
is desired these approximate methods should be abandoned and a fully developed 
mathematical program, computer solved, should be employed. 
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III. EXAMPLE OF PLANETARY SWING-BY FLIGHT 


The example presented here will be concerned with the flight of a spacecraft 
from an Earth orbit to (say) an outer planetary orbit. During this maneuver the 
spacecraft will pass close by an intermediate (massive) planet, executing a 
"swing-by" maneuver; thereby having its ballistic trajectory altered by the 
gravitational influence of this intermediate body. 

In this study the ballistic path of the spacecraft is obtained from a restricted 
two-body solution, with the sun serving as a primary mass. For this problem, 
then, the intermediate planet effectively acts as a perturbing influence, altering 
the two-body path, and 1 ’bending” the ballistic space trajectory according to the 
mode of passage by the disturbing body. 

In the analysis of this flight problem, two methods of approach will be pre- 
sented. For the first, and simplest, case the effect on the trajectory produced 
by the perturbing planet will be assumed to occur at a given, fixed radial distance 
(namely, at that heliocentric distance locating the disturbing planet itself). Thus, 
the spacecraft will fly to the perturbing planet's orbit where the ”swing-by” will 
be treated as a local occurrence — "bending” the flight path at this position in 
space - then, the spacecraft will be assumed to move on along a new ballistic 
trajectory. This is a rather uncomplicated approach to the "swing-by”; one which 
is easy to follow and to understand, yet one which provides the basic information 
needed to describe this type of a perturbed flight track. 

In the second case study, a more rigorous and complicated approach is 
undertaken. For this investigation a finite size of the perturbing planet's "sphere 
of influence” is assumed, and the "bending" of the space track is accounted for in 
a more realistic fashion. Here, the geometry of the (close-in) planetocentric 
space trajectory is determined, and evolves, as a finite dimensioned situation. 

In this regard the events arising from the perturbation are relegated to the 
"vicinity” of the planet rather than to a point in heliocentric space. The conse- 
quence of this latter approach is that the problem becomes more complicated, 
and needs additional analysis and computational effort; but it does lead to a more 
realistic evaluation of the physical problem. Whether or not the added effort, 
involved here, is worth the additional accuracy achieved is a matter to be decided 
by the investigator. His need for computational accuracy and for geometric 
description will dictate which of these two approaches (if either) should be 
undertaken. 
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III. A. AN INTERPLANETARY TRANSFER 


In this example it is proposed that an earth launched spacecraft wii 1 travel 
to one of the outer planets, aided (enroute) by a planetary "swing-by" maneuvei. 
For the sake of specifics, suppose that a mission is designed whereby Pluto is 
reached via a Jupiter swing-by, following the earth orbit launch. 

Having the need for an initial and comparison trajectory, then a first con- 
sideration will be to describe a basic ballistic transfer path; subsequently the 
swingby operation will be viewed as a modification to this reference trajectory. 

For a fundamental ref e rence path suppose that a Hohmann* transfer from 
(an assumed) earth circular orbit to an assumed circular outer planetary orbit 
is selected. Thus this basic flight path will be a nominal, least energy transfer 
arc, and the subsequent swing-by will result in a modification of this conic. 


III.B.l. THE HOHMANN TRANSFER 

As a means for describing the Hohmann trajectory, its several parameters 
will be listed balow — as these relate to this problem. It should be recalled that 
this trajectory is part of an ellipse which is cotangential to both of the assumed 
circular orbits (e.g. those of earth and the terminal planet). Thus, here, the 
earth's orbit radius will match the Hohmann pericentric distance while the size 
of the terminal circular orbit will correspond to the apocentric radius. 

The orbital (two-body/ >ar a meters: 

(a) eccentricity of the Hohmann ellipse: 


a — r + r 

peri rpo 


r peri + f apo 


r f “ f ® 
+ r f 


(b) major axis of the ellipse: 


2 an = r apo + r peri = F f + r ® 5 


<m.i) 


(in .2) 


*The Hohmann ellipse was selected simply as a matter of convenience. Any convenient initiol 
path could be chosen; because of its relative simplicity, in definition, the Hohmann trmsfer 
serves as on easily understood reference ore. 
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(c) focal parameter for the trajectory: 



r peri *"apo 


r . + r 

pen apo 



r 



+ r , 


(m.3) 


III.B.2. DYNAMICAL PROPERTIES OF THE HOHMANN ELLIPSE 
(a) the specific angular momentum: 


h H = | r x V 


r V 

pen peri 


* 


(measured at the pericentric position on the ellipse; the earth’s radial distance 
from the sun, 1 a.u.) wherein, 


F peri (^ r ©> = 


A Ph 


1 + e. 


and, from the specific energy equation, the speed at pericenter is 


V . = 

per i 



^ 

f peri ^ 



2r f /r ( 


© 


1 + r, /i 


1/2 


Here ^ is the gravitational parameter (Gm e ) for the sun. One should keep in 
mind the fact that this particular ballistic path is a heliocentric trajectory. 

Combining these definitions, it is evident that the momentum constant for 
the Hohmann path can be expressed as ; 


1/2 

(III. 4) 


wherein (circ ular or bit speed at a distance equivalent to the earth’s helio- 
centric radius) = /m q A e . 


= r e V c a 


2r. 


r f + r 6 
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(b) The level of energy, for this heliocentric trajectory, can be described 
from the specific energy expression; thus, 


g a / 

H \ ^ r /H ^ a H r © + r f 

or, after manipulation, 

-V c 2 

e. = — • (in.5) 

1 + *■, r. 

These characteristics, above, are adequate for describing the Hohmann arc. 
There is one parameter which should also be defined - the time for motion - this 
is deferred for the moment, it appears in a following section. 


III.B.3. CHARACTERISTIC SPEEDS AND VELOCITY IMPULSE SCHEDULE 

On the Hohmann transfer ellipse the two speeds which are of primary im- 
portance, immediately, are those corresponding to departure (at pericenter) and 
to flight terminus (apocenter, at the terminal planet radius). Manipulating the 
definition for specific energy it is easy to show these two characteristic speeds 
are: 


(\eri) H - (X®)h “ \ a 


2r. 


f f + 


1/2 


and 


( V apo) H = < V f) H = V_ 


2r a 


1 1/2 


f f + r € 


(in. 6a) 


(m.6b) 


wherein V Cf = /m© / r f . 

Next, in order to begin the transfer, at pericenter a AV (impulse) would be 
needed to shift the vehicle from its earth circular orbit onto the transfer ellipse; 
and, in a like fashion, to get off the eUipse (at apocenter) and onto a terminal 


circular path, a second AV is needed. In this regard, then, the AV schedule needed 
to achieve these two actions can be determined as: the impulse at, 


(1) departure from pericenter; 



and 

(2) terminus, at apocenter; 


(III. 7a) 


AV 


f 




(III. 7b) 


It should be evident that these AV. values are in the (±) direction relative to the 
local velocity vectors, respectively. 


III.B.4. THE TIME FOR TRANSFER 

The Hohmann transfer, which occurs over one-half of the ellipse, would re- 
quire a transit time equal to half of the orbital period; thus, the transfer time 
along this idealized flight path would be: 



or, using the Hohmann parameters, 


At H 



(III. 8) 


In the present problem this quantity is not needed since the full extent of the 
Hohmann geometry is not utilized, per se; however, it is included here since it 
is one of the parameters describing the reference path. 
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III.C. SPHERE OF INFLUENCE 


For those calculations which describe the "swing-by," per se, it is necessary 
to define the perturbing planet's "sphere-of-influence." This is the region in 
which the maneuver occurs. Consequently a ficticious boundary which defines 
the "sphere" can be obtained from: 


r 


IN 


= R 




(III. 9) 


wherein r IN describes the planetocentric radius of the "sphere of influence;" 
while R is the radius to the planet (m d ), measured from m 0 (the sun). Also, the 
quantities m 0 and m d are the solar and disturbing planet masses, respectively. 


m.D.l. THE SWING-BY MANEUVER; MODE I 

The simpler mode of analysis, used to describe the swing-by maneuver, will 
be designated herein as Mode I. For this operation the flight begins as (say) a 
Hohmann transfer, and continues as such until the spacecraft is in the very im- 
mediate vicinity of the "swing-by" planet. Then, in the neighborhood of this 
perturbing body - actually, within its sphere of influence - the swing-by occurs. 
Due to the planet's gravitational attraction the flight path is altered so that the 
subsequent heliocentric trajectory is no longer a continuation of the Hohmann 
ellipse, in general, but becomes another (and possibly different) heliocentric 
conic. 

In this first analysis the effects of the swing-by are assumed to occur at a 
fixed radius; that one which locates the perturbing planet in heliocentric space. 
Because of this assumption the swing-by is assumed to alter the trajectory 
locally, producing a change in the vehicle's heliocentric kinetic energy. However, 
the idea of this assumption does not allow any apparent change in the vehicle's 
potential energy (due to the encounter) . “ “““ ~ 

Now then, of necessity, when one does attempt to ascertain the change in 
direction for the spacecraft's velocity vector, an analysis which determines the 
flight path through the sphere of influence must be undertaken. 

In addition, the sphere of influence may be utilized as an aid in accurately 
estimating the transit time through this region of (planetocentric) space during 
the encounter. This particular topic will be discussed later! 
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III.D.2. A DESCRIPTION OF INTERCEPT WITH THE SWING-BY 
PLANET'S ORBIT 

The assumptions used for a Mode I swing-by consider the disturbance to, 
occur at the planet's heliocentric radius rather than in its sphere of influence. 
Consequently this intersecting of the paths of motion is an important factor in 
the analysis; and one which necessitates an evaluation of this intercept. In this 
regard one needs to know conditions, parameters, and a description of the 
spacecraft's path, at the swing-by planet's heliocentric radius (F d ). 


s'* PATH FOR md 

/ 



Figure III. 1-Sketch Showing Intercept of the Hohmann 
Path with the Orbit of m d . 

Since the characteristics of the Hohmann path are known, explicitly, then at 
encounter with the disturbing planet's orbit, the spacecraft's heliocentric position 
(angle) can be acquired from the conic equation. A manipulation of this expression, 
expressedly written to describe the Hohmann ellipse (at r d ) , gives 


(<P H )d = cos _1 



o - «i> - f • 

H _ 


(III.10) 


To define the spacecraft's speed, at ”m d ", use will be made of the specific 
energy Eqn. (solved by v H ); that is, the speed at r d , for the vehicle on its 
Hohmann track is; 
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(III. 11) 


(Vd = 


12 

i 







Also, for subsequent use, the heliocentric speed of the planet (m d ), at encounter, 
can be assumed to be 


V 


d 



(III. 12) 


since the path of M m d ” is presumedly circular. 

In order to describe the elevation angle (y Hd ), at r d , for the spacecraft’s 
velocity vector one could proceed in the following manner: 

In view of the conservation of (two-body) moment of momentum, write, 

V |K| 

V co *y«t = — = v 

d d r d 

Recalling that ^ = /p H ^ 0 ; then, accounting for the conic equation, rewrite 
the above expression as 


V H cosy„ 

H d H d 




(1 


+ € h COS CP H 



and, therefore, 


\ = c ° s '‘ 



(III. 13) 


(Note, there is no need to find the sin" 1 function (here) since the physics of this 
problem locates 7a d in the first or second quadrants of the orbit’s geometry. 
Similarly, any other transfer situation would physically describe y so that no 
ambiguity would exist.) 
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In view of the assumptions made for the Mode I solution the spacecraft meets 
the swing-by planet - and subsequently departs from it - at the heliocentric 
position (r d , cp H<J ). Notwithstanding this simplification, in order to evaluate the 
influence of the planet, m d , on the spacecraft, during the ’’encounter,” the overall 
analysis must include the essentials of a swing-by maneuver in order to ascertain 
the degree of trajectory perturbation introduced by this occurrence. For the 
purpose of determining a spacecraft track, relative to the encountered planet, 
the swing-by problem is treated (here) as a two-body situation, but with ”m d ” 
serving as the primary mass, now. 


III.D.3. TIME OF FLIGHT TO THE SPHERE OF INFLUENCE (MODE I) 

Before continuing with this analysis, and prior to considering the flight path 
inside the sphere of influence, it would be advisable to ascertain the time required 
to perform the first segment of the overall mission. This will account for the time 
lapse from the initial point, on the Hohmann arc, to the point of intercept with the 
m d -orbit (the perturbing planet's heliocentric orbit). It should be evident that the 
arc in question is from pericenter (at one earth heliocentric radius) to the inter- 
cept, at one m J heliocentric radius. 

d 

Knowing the position of this intercept (r d , cp Hd )» see Eqn. (III. 10); and knowing 
the orbital parameters for the (Hohmann) arc, then it is quite apparent that the 
time required to reach the intercept position is given by 


At i r - V, - t Q=J—r € h tan v + 2 tarrl 


1 €h . tan ^ W (III. 14) 


with e H ,a H described in Eqns. (.1, .2), and y„ obtained from Eqn. (III. 13). 


III.D.4. CONDITIONS AT ENTRY INTO THE SPHERE OF INFLUENCE 


At that point in time and space when the spacecraft reaches the ’’boundary” 
for the sphere of influence, it will be necessary to describe the trajectory con- 
ditions relative to the attracting mass (m d ). This is needed so that one can shift 
from a heliocentric track to a planetocentric flight path. 

At the imaginary boundary enclosing the sphere of influence, the vehicle's 
velocity vector(s) are related by 
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HOHMANN PATH 



Figure 1 1 1.2-Geometric and Kinematic Conditions at Intercept 
of the Hohmann Path with the Orbit of m d 


where V 2 is the velocity of the vehicle relative to m d ; V , is the velocity of m d 
relative to the sun (©); and, V Hd is the velocity of the vehicle, at the distance 
r d , on the Hohmann path. 

Solving for | V 2 | one obtains a measure of the speed for the vehicle, relative 
to m d , when it reaches the "sphere of influence." Thus, 

IV 2 l = /v^ d -2Vv Hd ; 

or, noting the previous sketch, 


V 2 = V V d + -2V d V H cosy H . (III. 15) 


III.D.5. THE TRAJECTORY INSIDE THE SPHERE OF INFLUENCE 

Once the vehicle enters the sphere of influence its flight path is considered 
to be a planetocentric two-body trajectory. It is most likely, here, that the pre- 
viously described heliocentric ellipse becomes a planetocentric hyperbola. Of 
course this will depend on the magnitude of V 2 in relation to (say) the planeto- 
centric circular speed (or escape speed), at the radius of entry. 
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Inside the sphere of influence, where the only gravitational attraction is 
that due to mass m d , the flight path is (again) characterized by constant specific 
energy and moment of momentum. In view of the assumptions used for the Mode I 
calculations the planetoeentric study will treat this part of the problem as if the 
entry conditions were the same as the "infinity" values for a two-body hyperbola. 
That is, from the energy expression in planetoeentric space, 


V 2 

T’T " 2a 


(which presumes a hyperbolic geometric figure), the energy relation - correspond- 
ing to large radii - leads to the approximation 


2 ” 2a 

Allowing this approximation to be a description of entry into the "sphere of 
influence ," it follows that the planetoeentric trajectory has, as its major dimen- 
sion (approximately), 


a 



(III. 16) 


vith V 2 defined by Eqn. (III. 15) above. 

As an aid to obtaining pertinent information for this "hyperbolic encounter," 
vith the perturbing mass (m d ), one can employ (first) the conic equation. Thus, 
i description of the planetoeentric position angle, relative to the planetoeentric 
>ericenter, at entry into the sphere of influence, will come from 


Pd 

r 

1 + e d cos <p tt 

'here qp* implies the "infinity," or large radius, (entry) position. Following with 
le idea that r x is very large, then 


~ 0 = 1 + e d cos cpgjj 
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(hence cp m essentially defines the approach - and/or receding - asymptote); 
or, in terms of the complimentary angle (w), where \y = v - qp^ (see the sketch 
below), 





Figure III. 3-Slcetch Illustrating the Planetocentric Hyperbola 


JL - 0 = 1 + cos (rr - y) , 

and 


1 

cos vy = + — • 
€ d 


(III. 17) 


From this expression the eccentricity of the planetocentric conic can be obtained 
as: .. 


e d = sec \y • (III. 18) 

A geometric parameter of major importance to the swing-by is the angle 'V' 
shown on the accompanying sketch. This angle describes the magnitude of the 



rotation which is given to the spacecraft's velocity vector as a consequence of 
the "swing-by” maneuver. Denoting the vehicle's final velocity vector as V 3 in 
planetocentric space - just before it departs from the sphere of influence - then 
the angle between V 2 and v 3 is "a”. Thus, in mathematical parlence, 


a - cos -1 



(III. 19) 


also, it is the angle which describes the "size" of the "kink" produced in the 
heliocentric space trajectory due to the "swing-by" operation. 

A description for a can be acquired in the following manner: Since a = 
<p w - w (see the sketch on preceding page), then 

cos a: cos(cp tD - \|/) = cos cos + s in cp^, s in 'V ; 

and 


sin a = sin cp^ cos \p - cos sin 
hence, after substitution and clearing, 

€ d " 2 2 / e d - 1 

cos a = j — , sin a = , 

e d £ 5 

taking into account Eqn. (III. 17); and/or 


(in. 20a) 


l/e\ - 1 

tan a = d — . (III.20b) 

4- 2 

A study of this last expression will indicate that the angle a is rather strongly 
dependent on the eccentricity of the planetocentric path - and, as an a priori 
statement, the figure is a hyperbola. Also, it is apparent that t d , V and a are 
related according to the tabulations noted on following page: 
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* h|1 JJpL^ r " 





When 

€ d = 

Values for 

y = 

Values for 

U SS 

1.0 

0 

7 7 

1 < e d < 

0 < \\i < 77/4 

77 > a > 77/2 

A 

77/4 

77/2 

> /2 

77/4 < < 77/2 

77/2 > a > 0 


The "size" of this hyperbolic path about the disturbing mass (m d ) can be de- 
fined by means of its focal parameter (p d ). That is t since 

vj- »(«J - 1); 

then according to the a ssumptions set forth for this pro blem’s study - namely, 
that V 2 ^ = /f± d /a\ also, since tan a/2 = - it follows that 


P d - 2 (HI.21) 

(v, *«■ fT 

In general, the characteristics of the planetocentnc (hyperbolic) trajectory 
are known, now, and this geometry is described, at least in principle. 

The conditions (speed, orientation of V , etc.) at exit from the sphere of 
influence have not been determined but will be obtained below. Also, even though 
the size of the sphere of influence is not considered - on the heliocentric scale 
of distances - it may be desirable to estimate the time required to move through 
this figurative region in space. To this end, in a subsequent paragraph a "time 
needed for the planetocentric transit" will be defined. 


HI.D.6. FLIGHT CONDITIONS AT EXIT FROM TH3 SPHERE OF INFLUENCE 

The assumptions made for the Mode I swing-by computational procedures 
dictated that the discontinuity in the spacecraft's trajectory, produced by the 
encounter, would occur at the fixed position locating the disturbing planet in 


.7 



Ipsi^p 


heliocentric space. This would imply a vanishingly small sphere of influence 
(on a heliocentric scale); however, in order to view and describe flight conditions, 
at exit from the sphere, it is necessary to give this region some finite size (as a 
planetocentric space). In this manner the orientation, etc. of the vehicle's "exiting" 
velocity vector can be defined. The pair of sketches below will indicate condi- 
tions which are of interest here. 


i 




VELOCITY DIAGRAMS 




HYPERBOLA (FROM • TO m<j) 


Figure III. 4-Geometrics for the Encounter at m d 


The sketches (above) illustrate geometric conditions which fit the ideas of the 
Mode I calculations in the region of the encounter with m d . For the sake of clarity 
the nomenclature used here is described below: 

(1) At the position of m d , the heliocentric velocity for the vehicle is V Hd ; the 
velocity of m d (itself) in this space is v d ; also, according to the following vector 
expression these are related by 


with V 2 being the planetocentric velocity of the spacecraft at entry (A) into the 
Sphere of Influence. 
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(2) Point (B), at exit from the Sphere of Influence, is noted to have the 
following vector relation satisfied: 

v 3 + v d =v Td . (in. 22) 

Here V 3 is the planetocentric velocity of the spacecraft, and ^T d is the helio- 
centric vector (the subscript T refers to the terminal heliocentric arc ) . 

(3) In heliocentric space the elevation angles (y H , 7i d ) locate V Hd and ^r d 
relative to V d (the planet's local velocity vector, paralleling the heliocentric 
horizon) . The corresponding planetocentric angles are (y Hd + S Hd ) and (y Xd + S Td ), 
locating the planetocentric velocity vectors (V 2 , V 3 ) relativeto V d . (One should 
note that the S. k angles position the V 2 , V 3 relative to the V. k vectors). 

(4) Since the hyperbola is a symmetric figure in planetocentric space, then 
it follows that 



(III .23) 


Also, according to the assumptions in use here the heliocentric speed of m d 
is described by 


V 


d 



(5) In planetocentric space the vehicle flies a hyperbolic track, about m d , so 
that the effect of the encounter is to alter the direction of the spacecraft's velocity 
vector by an (angular) amount, a - in that space. 

In heliocentric space the influence of the encounter is to alter the direction 
of V (spacecraft velocity) according to the following scheme: Relative to V d the 
approach velocity (V„ d ) is located by the elevation angle ( y Hd ) . After the encounter 
the angular position of the vehicle's heliocentric velocity vector (V Td ), relative 
to V d , is (y Td ). 

From the sketch (above) it is evident that V 2 is described relative to V d by 
two relations; namely, 

* <V V d ) = r„ d + 8„ d = y Td + S Td + a. (Ill ,24a) 
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Next, in heliocentric space the angle between the vectors *V T and V H ^ is 

* ( \ ’ V = * 7 »d " 7l d ' ; 
or, from the expressions above, 

* (V Td , V„ d ) = | a + S ?d - S„ d 1 , (III. 24b) 

denoting that the velocity change in direction is not the same in heliocentric space 
as it is in planetocentric space. 


III.D.7. A DESCRIPTION OF THE VELOCITY, AFTER ENCOUNTER 

At exit from the sphere of influence the spacecraft has a velocity V 3 , relative 
to m d , and a velocity V Td in heliocentric space. (Note: For clarity the subscript 
”T” will be employed hereafter to identify the post-encounter, or terminal, 
heliocentric arc.) 

This heliocentric vector can be fully described by: (1) defining its magnitude; 
and, (2) indicating its direction (say) in terms of its elevation angle (y Td )» The 
speed (V Td ) is acquired by manipulating Eqn. (III.22); that is, 

v Td £ |v T J = + Vj + 2V 3 • v„ , 

or 

\ = y n + 2V 3 V d COS (y Td + S Td ) , (III . 2 5 ) 

wherein (as noted previously) 


V 3 VT-- a " d V d = 



also, the angle (y Td + S Td ) is described from Eqn. (III. 24) in terms of a, y H , S H . 


In order to describe y Td , explicitly, the geometry shown in the sketch above 
should be consulted to ascertain the following relations: 
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On the "velocity diagram" one notes that 


V3 s i n ( § x + y T ) — Vy sin ( ) 


but (8_ +%.)*[ (8 + y_ ) - a] from Eqn. (III.24); so, 

T d T d H d T d 


V sin [(5 +7 ) -a] 

sin (y T ) = 2 2 


(III. 2 6a) 


and, by manipulating Eqn. (111.22), one can show that 


(V? +vj)-vj 

cos(y T ) = — 2 

d 2V_ V, 

T d d 


(III. 26b) 


A combination of these two expressions leads directly to the arc tangent expression; 
or, 



2 V d V 3 sin f(^H d + 




(111.26c) 


Equations (III.25) and (III.26) should provide sufficient information to properly 
define the velocity (V Xd ) in heliocentric space. With this as a known quantity, 
plus a knowledge of the position ( 7 d ), then according to previous assumptions one 
should be able to establish the nature of the terminal arc. Having a description 
of this (last) segment of the space track, the analyst would be able to complete 
his computations for the entire mission. 


Prior to beginning a description for this last phase of the journey the side 
issue of "time to fly through the sphere of influence" will be mentioned. 


m.D.8. TIME OF FLIGHT THROUGH THE SPHERE OF INFLUENCE 

In many instances it may be desirable to estimate a time required for the 
spacecraft to complete the hyperbolic encounter, per se. Recognizing that this 
is the time of transit through the sphere of influence (from point (A) to point (B) 
shown on the sketch above) , then it is readily evident that the time lapse for this 
segment of the flight track can be defined from 
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(III. 2 7) 



c d V f d - 1 sin ^2 

1 + 6 d COS Cp 2 



where (At) IN signifies the time lapse in traversing the Sphere of Influence. In 
this expression (e d , a) are parameters describing the hyperbolic arc; M d is the 
gravitational parameter for the perturbing planet (m d ); and r P 2 is the planetocentric 
position angle corresponding to the actual entry (position) into the sphere of 
influence. The multiplier, ”2”, is included here to describe the full transit time ; 
Eqn. (III.27), without the multiplier, is indicative of the time of motion from 
pericenter to the sphere of influence boundary (described by qp,). 


In order to obtain a realistic time measure, here, it is necessary to relax 
the assumption-used earlier- regarding the "size” of the sphere of influence. 

In the present instance, a true value for r IN is needed so that a proper estimate 
of <p 2 (the position angle, corresponding to r IN , on the hyperbolic path) can be 
obtained. Thus, using the conic equation - naving determined r IN from (ITT. 9) 
and e d from Eqn. (III. 18) - it follows that 


<P 2 = 





with p d found from Eqn. (III.21). 


(III. 2 8) 


In keeping with the various assumptions used for the Mode I calculations, it 
is quite likely that for the simpler case studies, the time to fly through the sphere 
of influence could be neglected (in many instances); yet, on a small (time) scaled 
mission it may be prudent to include this time fraction in the overall estimates. 
Such a decision must be made as a part of an evaluation for the overall problem. 


III.D.9. CHARACTERISTICS OF THE TERMINAL TRAJECTORY 

After the spacecraft exits from the influence of m d , and continues on its 
way to m f (the "final mass" and its space position), via a free flight heliocentric 
trajectory, one should ascertain the characteristics of this flight path. As an aid 
to this description the flight path can be identified simply in terms of how the 
vehicle’s heliocentric speed, after swing-by, compares to (say) the heliocentric 
escape speed. If the spacecraft’s speed is less than the escape value then the 
final path is an ellipse. However, if the flight speed is above the parabolic value 
(corresponding to the local value of r ) then the vehicle escapes from m d and 
progresses on toward m f along a hyperbolic "free path." 
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In order to describe the character of this final path the vehicle's speed con- 
dition may be examined according to (the specific energy equation): 
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but/x /a = V 2 , so that 
1 © 


V 2 r V 2 _ v 2 , (ni.29) 

wherein V| #Co = 2 /z 0 /r d (the escape speed squared, in heliocentric value), and 
V a>o is the heliocentric hyperbolic excess speed. 

If Eqn. (m.29) is negative then the "escape" condition has not been realized 
and the flight path "T" is elliptic; however, if V 2 >0, then escape has been 
attained and the heliocentric trajectory is a hyperbola - incidentally, one should 
recall that the parabolic path is a special case of the hyperbola, but one where 
a T - oo and V* vanishes. 

Necessarily, the Terminal Flight Path would not be expected to have the same 
(or a parallel) line of apsides compared to the Hohmann Path. Likewise, these 
figures (in heliocentric space) are not likely to have a same eccentricity, or 
size. So, in order to complete the description of the terminal path these orbital 
parameters must be obtained. 

The eccentricity of the terminal path (e T ) is easily shown to be expressed by 


e 2 = 1 + 


2Eh 2 


or 


e T - 


r V 2 

a V T 


1/2 


— - 1 \ cos 2 7 t + sin 2 y_ 
*d r d 


(III. 30) 


for the situation considered here (Mode I transfer geometry); and at exit from the 
sphere of influence, but with the exit assumed to occur at r d . 

The size of the terminal path can be described by the parameter (p T ); or, 
by the major dimension (a T ). Since |E l | = /x /2a (for either figure), then it 
follows that 
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ill 

© 

=L 


2E, | 

VJ - 2 — 

T d r d 


(III. 31) 


where the absolute value is employed to leave a T as a positive definite parameter 
regardless of the geometry describing the path (T). 

In order to account for the ’’shift” (angular displacement) of the line of 
apsides (between the original Hohmann, and the terminal path) one should calculate 
the (new) position angle (cp T ) and then describe this displacement according 
to the relation d 


A<p d = <p„ d - cp^ , 

where <p h is obtained from (III. 10), and where cp T is obtained from 

d d 


(III .32) 


<P T = c 

T d 


os 1 r pt _ i . 

L r d e T e T 


(III. 33) 


In Eqn. (III.33) p T - a T | 1 - e* | , accounting for all conic types (except the 
parabola). 

Note: If (III.32) yields a negative resultant, then the line of apsides (terminal 
path) has moved in a retrograde direction relative to the original apsidal axis 
(for the Hohmann). A positive resultant indicates a posigrade rotation of the 
apsidal axis (wrt the Hohmann). 

Some other quantities of interest, describing the terminal path, are noted 
bolow: 

(1) Pericentric radius (of pathT), 


r > T = a i I 1 - S' 5 


(2) Apocentric radius, 

v T r ®i o + £ t> 

(3) Period of Motion, for path T , 


(hyp. ellip.) 


(ellip. only) 


p ■ : 


(ellip. only) 
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(4) parameter, p T , 


P T - a T 



and 


(5) parameter p T (parabolic path), 

Pt * 2, > t • 


(elltp., hyp.) 


(parab. only) 


| 

f 

I 


m f ORBIT 



\ 

Figure III. 5- Sketch Showing the Flight Paths for the Swing-By at 


1 


At the terminus of the mission the spacecraft is assumed to have reached 
m f (along the Terminal Path). In this regard the position of the intercept at the 
terminus is specified by (r f , <p f ) where r f is known, a priori. To calculate 
q>_ one can resort to the conic equation, obtaining 

T f 


<p_ = cos 1 

T f 




(in .34) 


wherein p T and e T are known (see Eqns. (HI.30, (ni.33)). Note that Eqn. (m.34) 
locates m f relative to the pericenter for path (T ) . 
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III. D. 10. TIME OF FLIGHT SUMMARY 


In the analysis discussed here the total time of flight involves motion along 
two arcs , at least. The first segment involves the motion from © to m d ; and, the 
second involves the arc from m d to m f . (In some instances, but not necessarily 
the case at hand, one may wish to account for the time within the sphere of 
influence . However, for present purposes this latter time increment will be 
ignored, and the mission, timewise, will be assumed to be composed solely of the 
two principal arcs for the transfer.) 

The two arcs considered here may be segments of ellipses; or, the last arc 
may be a segment of a hyperbola. Thus, to account for the total ’’time of motion” 
it will be necessary to include both (possible) equations and to use them appropriately. 

Since the equations for time, set down below, count time from pericenter 
passage it will be necessary to ’’manipulate” these expressions in order to obtain 
the desired results. One should note that these expressions involve the orbit 
parameters (a, e ) and the position coordinate (cp), only; thus, time is expressed as 

t = t (a, e, op); 


or specifically: 

(1) For an elliptic arc (from pericenter to some position cp ); 


t = 



J e 2 - lj sin cp 


1 + e cos cp 


+ 2 tan 


Vt 


- e cp 

— — t an 

+ € 


(III. 35 a) 


(2) For a hyperbolic arc (from pericenter to a position cp ); 


t = 





1+6 


lsin<p .-2tanh-» (111 I 

COS Cp \ ? € + 1 


tan I 
2 


(III. 35b) 


or 




v/e 2 - 1 sin cp 


- t n 


1+6 COS Cp 


tan(cp/2) 
t an(cp/2) 


. (III.35c) 
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In computing the time of flight (neglecting the time passage when the vehicle 
inside the sphere of influence) the above equations are employed in the following 
manner: 

(A) Time from © to m d - use (III.35a) wherein 

t = t (ajj, <p h , € h ) - see eqn. (III. 14), also. 

d 

(B) Time from m d to m f - use (III.35a), or (m.35b), (III.35c), as follows: 

(1) if 0 < e T < 1.0; then use (III.35a), calculating: 

tj = tj (Sp, 

t 2 = t 2 (Sp, CP T , e T ), 
d 

and then, 

At = tj - t 2 = time of flight from m d to m f ; (III. 36a) 

(2) if 1.0 < e < oo ; then use (III. 35b) or (III.35c), calculating: 

*1 = l l C 0 !’ %■ € T>’ 
t 2 = *2 ( a T € T^ 

and then, 

At = tj - t 2 = time of flight from m d to m f . (III. 36b) 


III.D.ll. ENERGY CHANGE DUE TO SWING- BY 

Another of the trajectory parameters which can be significantly affected by 
the swing-by maneuver is t v j total specific energy, referred to heliocentric space. 
The change which occurs in the energy is caused by the vehicle's heliocentric 
velocity vector being altered (in direction), due to its passage closeby the dis- 
turbing body (which is, itielf, in motion). 

In order to describe the change in total energy (positive or negative) one should 
recall that the specific energy is defined by 
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''iMwWlll,...,,; 


thus, in heliocentric space, V is the speed along the arc(s) (H and/or T), while 
/i ~ /x @ , and r defines position relative to m e . 

For this present case, the position at encounter has been assumed to be fixed; 
i.e. r , and consequently, 



or 

AE = Jl - J!i . (IK-3 7 ) 

2 2 

This last expression will describe the energy gained (AE > 0), or lost (AE < 0), 
during the swing-by maneuver. 

Essentially, now, the calculation procedures, by a Mode I analysis, are com- 
plete. In the next paragraphs, the more exacting scheme (Mode II) will be described. 

III.E.l. MODE II TRANSFER, SWING-BY OPERATION 

The extension to the earlier problem, introduced here, is such that a finite 
size for the sphere of influence is included in this next formulation, and is 
utilized throughout the computational procedure. That is, instead of having the 
two principal arcs end and begin at m d , these are now carried only to the sphere 
of influence - where each arc intersects this imaginary boundary, and the motion 
through the sphere of influence becomes an integral part of the problem. The 
motion through this region becomes a real and essential piece of the overall 
analysis. 

The assumption of r IN as a finite (or non-zero) value adds a complication 
to che problem which was not previously encountered. 

In describing the sphere of influence, its size will be assumed to be calcu- 
lated from (see Eqn. (III. 9)) 


P i 


r. 

1 



» 
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where i , j refer to the disturbing mass bodies being considered; or, specifically, 
for this problem where the encounter is with m d , then 


/5 

, p l - radius of sphere of influence, (III. 38) 
relative to m d as a center; 

and, wherein r d locates m d relative to m 0 . 

Now, the Mode II problem will be concerned with three (3) separate arcs of 
motion, namely; 

(1) The Hohmann arc from ® to (p,) H , the sphere of influence; 

(2) the hyperbolic arc through the Sphere of Influence; ami, 

(3) the free trajectory from the Sphere of Influence to m { ! 


m.E.2. THE PRE-ENCOUNTER TRAJECTORY 

This path arc is (again) selected to be a portion of the Hohmann trajectory* 
connecting the orbits of and m f . As expressed in Eqns. (III.l) through (III. 8) 
the characteristics of this arc are: 



( a ) € h - 


^■po ** ^peri 
*"apo ^ *"peri 


r f “ r ® 


r f + r ® 


(eccentricity of the Hohmann arc) 


0>> V H 





(speed, at any point along the Hohmann arc) 


(c) 




(speed, at position on 
the Hohmann arc) 


*As noted earlier, the selection of o Hohmann path is done mainly for convenience and familarity. 

In any case the initial arc moy be selected as any conic the analyst would care to describe. 
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<d) 





* r< 


© 


(speed, at the m f position 
on the Hohmann arc) 


(e) AV. V 
(t) AV f --V Cf 



(increment in speed needed 
to leave the circle ( rj and 
get onto the Hohmann arc) 

(speed increment needed to 
get off the Hohmann arc, at 
m f , and onto the circle ( r f )) 


(g) P H a H (1 f 

(orbital parameter) 

( h ) b H '■ *h Z 1 “ c i = V^Ph 

(minor axis length) 



(time for motion, from m. 

w 

to m f , on the Hohmann 
arc) 


The quantities listed above describe the Hohmann ellipse (or, specifically, 
the semi-ellipse) on which the spacecraft is assumed to move initially. What 
has not been described here is the location and conditions of motion at the inter- 
cept between the vehicle's path and the "boundary” for the sphere of influence. 


Because of the non-zero size for this boundary the "point" where intercept 
can occur is rather vague; actually, there is a fairly large length (of the major 
circle - this sphere) on which this point can lie. Assuming that the situation 
to fc described (uere) is fully one in a single plane, then the intercept can be 
described, and the trajectory characteristics at intercept can be determined. 


II1.E.3. SPEED AT ENTRY TO THE SPHERE OF INFLUENCE 


According to Figure III. 6, it is seen that the maximum magnitude for , 
r Imax = r d + /'j ; and, the least magnitude if r Imin = r d - p l . In a like manner, 
the range of values for the angle cp„ is; 

H i 



(III. 3 9) 
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Figure 111.6— Sketch Illustrating Extremes for 
the Sphere of Influence . 


where 


(these angles (cp. ) are measured relative to the pericenter of the Hohmann arc, 
with m 0 as a focus). 

For the situation illustrated on the figure, it is evident that 


s cp | > 0, S(p d < 0 ! 

•max min 

Note: At Sq) d | maXimin the radius, Tj , is tangent to the Sphere of Influence so that 
the intersection of the Hohmann trajectory with a major circle of the sphere is 
such that | "Fj | < 1 7 d |, by some "small” amount in a real sense. It should be 
evident that at this condition 7j and p. are orthogonal vectors (the scalar product 
vanishes); and, that, as such, 


r 


i 




9 


while, for this condition, 
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and the relationship between 6 and Sep is such that 


6 + I Sep 


max , 
min 


^ 77/2 • 


This would indicate that for this condition 


cos 9 - sin Sep 


( S ep = S ep 


max , mi n 


)• 


III.E.4. A DESCRIPTION OF THE INTERCEPT POSITION (I) 

The intercept to be considered here occurs as the vehicle flies along the 
initial transfer path. It is to occur in heliocentric space at the position denoted 
by <p and measured from the pericenter position for the Hohmann ellipse. 
Necessarily ^ satisfies the constraint, noted above; namely, that 

= <P Hd ± S<p d , (III. 40) 

where ± §cp d < sin (p l / r d ). Following along with this condition, then, it is 
apparent that the radius to the point, I , can be expressed as 


r 


I 


Ph 

1 + e H cos <p H 


(III. 4 la) 


from the conic equation. 

As an alternate description for this radius, note that from Figure III. 7 (for 
the condition shown) 


so that 


sin 8<p d = Pj sin 6, and, r r cos Scp d + p 1 cos & = r d 


Pj sin 6 r d - p l cos 0 

i-j = ; and/or, r x = 

sin S<P d cos 8cp d 


(III. 41b) 
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Geometric and Kinematic Conditions at Entry into the Sphere of 













In Eqn. (III.41) the radius (r x ) was defined; consequently, now, one knows the 
coordinates locating the point where the vehicle penetrates the ’’sphere of influence.” 
Even though this is a ficticious boundary, nevertheless it is essential that the 
intercept point (I) should be specified in order that the problem may proceed 
according to the assumptions set down earlier. 

Knowing the position of (I), one can describe the spacecraft's heliocentric 
speed, at this point, from the specific energy equation. In this regard, it is easily 
verified that 


d 


A 






1/2 



(See Fig. III. 7 for a graphical description of V H !) 

d 

Also, from the figure, it is evident that 


thus 


+ P T = r . 


(III.42) 


and, therefore, 


r i = r d - Pi ; 


r ? = r d + Pi ~ 2 r d • Pi 

= r d + P\ ~ 2 r d Pi cos ( e + 5 (ip d)- 
Including Eqn. (III.41) into the above relation, then, 

P H 

TJ = r d + p\ - 2 r d Pi cos (° + 8( ?d )• 

(1 + e H cos cP Hi ) 

wherein the unknowns (here) are the angles 6 and 6cp d ( r d is known, a priori; 
/Oj is calculated from Eqn. (III. 3 8); p H and e H are calculated for the Hohmann 
path; and the angle cp„ = cp u - Sep . is as described before). 

H i H d d 

A part of the present computational dilemma can be overcome by noting 
that 6 and 8 rPd are connected through Eqn. (III.41b) - which can be used to 
express 0 = A(8(p d ); that is, 
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i 


9 - sin' 1 


rj sin Scp d 


- sin 


-l 


Ph 

1 + e H cos cp H 


'sin Scp d 

~TT 


(III .43) 


Now, with 9 defined, as shown in Eqn. (III. 43), then 


Ph 


[1 + e H cos(cp H + 8<p d )] : 

d 


= + p\ - 2r A p l cos (6 + Scp d ), (III.44) 


which is an explicit relation in the unknown Scp d (i.e. including Eqn. (III.43) to 
relate 9 to S<p d ) . Of course, it should be remembered that Scp d is constrained 
as indicated by the comment given with Eqn. (III. 40). 

An iteration and/or interpolation formula (or algorithm) can be developed 
to satisfy this expression, Eqn. (III.44). It should be evident that there are 
actually four values for r which could theoretically satisfy this relation (based 
on a given value of | S<p d | ); these represent the various quadrant intersections 
of the path (H) with the sphere of influence boundary at a value of (±) Scp d . 
Necessarily, the physical problem and its constraints will dictate which of the 
points one should consider. 

After having obtained Scp d from (III.44) - and, likewise, defining 9 by Eqn. 
(III.43) - then the entry point is known, specifically and in toto, by the coordinates 
( r j, <p Hl ). (This is tantemount to specifying a value for q> H and then solving 
for r x ). 1 

A last quantity needed to fully describe conditions at intercept with the sphere 
of influence is the velocity vector’s elevation angle (y„ ) as noted on the figure. 

Mj 

By definition this quantity is obtained from 


e„ s i n cp„ 

tan y H = H Hl , (III. 45) 

d 1 + e H cos rp H ^ 

which is used to "locate” the velocity vector relative to the local horizon, (which 
is normal to here) . 

Occasionally one wishes to know the angularity of the velocity vector (V H ) 
with respect to the radius vector (p t ) - see Figure III. 8! From the figure it Is 
noted that (per the case shown) 
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Entry Conditions; Showing Planetocentric and Heliocentric Velocities, Angles of Reference and Coordinates 



* ( V 7 > ) = 7'V 


and, also, that 


77 


Pi) = (^+ l*dl>- ( y " >« d 


= ^+ I s^ d I +r Hd -yi 


(111.46) 


thus the angularity is defined in terms of known quantities. 
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* 
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III.E.5. TIME OF FLIGHT (FROM LAUNCH TO INTERCEPT, I). 


Since the initial segment of the flight path is an arc of the Hohmann ellipse 
then one can calculate the time required to go from pericenter (at m e ) to the entry 
point (I) by means of, 



(III .4 7) 


wherein tan y ttd is defined above, Eqn. (ni.45). This is only the first time incre- 
ment for the overall flight. The remaining two (2) intervals, to be calculated, 
will be discussed and determined subsequently. 


m.E.6. HYPERBOLIC TRAJECTORY, INSIDE THE SPHERE OF INFLUENCE 

Within the sphere of influence the flight trajectory will (usually) be a hyperbola 
thus the reference in the literature to a ’’hyperbolic encounter.” In this section 
the nature and description of this part of the overall flight maneuver will be pre- 
sented and discussed. 

On Figure (III. 8) the geometric and kinematic conditions at entry into the 
sphere of influence are depicted. The velocity vector for the vehicle, relative 
to the heliocentric (i.e. V Hd ) and relative to the attracting mass (i.e. V 2 ), are 
shown, along with the velocity of the disturbing mass (nid) relative to the helio- 
center (i.e. v d ) - and these are related to one another by 
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(Ill .4 8) 


V H d =W 


Accordingly, this expression can be scalar expanded to yield 


v 2 d = v 2 * V3 ♦ 2V a • v d = V* + V» + 2V 2 V d cos(V 2 , V„), 


(III.49) 


where 


i ( V 2 > v d ) - ( S H + >h + I S<P d I ), as shown, 
d d 

Rearranging Eqn. (III. 4 8) to read 



then as expression for | V 2 1 is obtained as 


(III. 50) 


V 


2 



V2- 


2 V, 


Hj 




(III. 51) 


wherein $ ( V H d » V d )-(7H d + I 8( P d |), which is known! Consequently Eqns. (III. 51) 
and (III.49) can be used to determine | V 2 | , in addition to the angle S Hd (= $ (V 2 , 


m.E.7. A MEASURE OF ASSUMPTION ACCURACY 

One of the assumptions used here is that the hyperbolic trajectory through 
the sphere of influence is reasonably close to a two-body conic and its properties. 
One recalls that for such a geometry the asymptotes and the figure merge at large 
distances from the focus (as the radius approaches infinity) . Certainly this 
situation does not exist here; and, even though the sphere of influ mce may be 
large there will always be some small (angular) separation between these two 
geometric lines. One measure of how accurate the present analysis is, compared 
to the assumption stated, would be the proximity of the radius (p x ) to the asymptote. 
That is, if the angle (shown below) vahished, then the approximation becomes an 
’’exact” hypothesis; and, as /3 becomes large (compared to zero) the degree of 
approximation is worsened. In this regard an expression for /3 is developed 
below; however, no numeric is assigned to the angle, as a means of ascertaining 
the degree of accuracy achieved by any particular analysis. 
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PO-a cos 4b *“7 
FOR/3 «o: 



Figure III. 9-Sketch Showing the Planet Centered Hyperbolic 
Path through the Sphere of Influence . 


According to the sketch (above) the perpendicular to the asymptote (drawn 
through the Focus F d ) - i.e., the line F d Q ( - b) , - is related to p (the radius 
of the sphere of influence) by 


F d Q(sb) = p 1 sin p 


(III. 52) 


where f3 = $ (V 2 , p x ). From an inspection of Figure III. 9 it is found that 

p = (V y H d + '^-(i - 6 ) 

= s t y + | S(P d | + 0- i , (III. 53) 

d d 2 

wherein all quantities are (now) known. 

An alternate evaluation of p could be achieved by determining the angle (w) 
of the asymptote, and the position angle (<p) corresponding to p v , then describing 
P accordingly. This should be compared with the expression above to ascertain 
the deviation between the exact and approximation situations. 
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III.E.8. DETERMINATION OF PERICENTER 


A quantity necessary to the successful study of an encounter problem (in a 
physical se.ise) is the expected planetocentric pericenter radius, pp. In order 
to relate this quantity to "entry” conditions, and known terms, one might follow 
the procedure outlined below: 

From conservation considerations (assuming, as noted before, the autonomy 
of the sphere of influence); 


and 


Pp V p 2 (F d Q) V 2 ^bV 2 , 



2 0 - = \ 2 - 2 

Pp 2 


Pi 


Combining these expressions, elimination Vp, then 


~ pp) ~ <% 


leading to a definition of pp as, 



III.E.9. ORBITAL PARAMETERS (HYPERBOLA) 

The orbital parameters for the hyperbolic path are easily defined from the 
geometric, dynamic properties of the encounter (within the sphere of influence). 
These are evaluation below: 

(1) The major length of the geometric figure (a) can be obtained from 
the planetocentric energy expression (E } ), since 



Md 

r 


+_— (assuming a hyperbola) 
2a 
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thus, for instance, for conditions corresponding to entry at the boundary 
of the sphere of influence, 


(III. 5 5) 


(2) Since the figure of the motion is a conic, then by definition, the radius 
to pericenter is 



Pp = a(c - 1); 


consequently, the eccentricity can be described from, 


e = — + 1 . 
a 


(III. 56a) 


As a matter of interest, comparing the expression above, for p ]d , to that 
given in Eqn. (in.54) it is noted that (alternately) 


e 



(III. 56b) 


(3) The turn (or angular change in direction) fcr the vehicle and its velocity 
vector, due to the encounter, is indicated on the figure above by a . This 
angle can be defined from the following: 

Note that 


a + 2 y = 77 , 


wherein the angle has been obtained from 

. b . 
y = tan 1 — ♦ 
a 

£ 

or, alternately, 

i 

I V ~ TT - 


1 

| 
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with 


. i 


m 


(recall that is the limit position angle for the hyperbola). If one employs 
these approximations, then 

a - v - [where = tan'^b/a)! ; (III. 57a) 


or, alternately 


o = - 


(III. 57b) 


This w'ould indicate that a y - / 3 for the encounter (regardless of the initial 
situation), ooviously this leads to an anomalistic condition which should be 
recognized as having arisen from the fact that, herein, the hyperbolic figure 
hap been assumed to lie close to its own asymptotes; and, that the entry 
conditions to the sphere of influence closely approximate the limit conditions 
for a (planetocentric) hyperbolic trajectory. Necessarily, the degree of 
error which has been incurred in this analysis is directly linked to the extent 
which the initial entry' state departs from the hyperbolic limit state. Some 
measure of this approximation can be linked to the angle, p • 


A check on che accuracy of w may be afforded (in part) by the following 
manipulations: 

(a) Knowing p. , then a value for is obtained from the conic equation 
(exactly) as, 


W im (=<P a ) = COS 


-1 


a(e 2 - 1) _ 1 
ep T ~ 


then according to Fig. III.10, 


* = ” - * ft' 

where p has been given by Eqn. (in.45)! (b) Using the result above, calculate 
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a = 7r - 2 W , 


and compare this to Eqn. (III. 57a) - or to (III. 57b)! 

With the information which has been obtained to this point in the analysis one 
is ready to proceed, next, to a description of the exit conditions from the sphere 
of influence. 

The flight through the sphere of influence has occured over a figure which 
has dynamic and geometric symmetry; thus the speed at entry (V 2 ) is the same as 
the speed at exit (V 3 ) - in agreement with the assumptions made earlier. That 
is, 

|V 3 | - |V 2 |. (III. 58) 

Also, in agreement with the assumptions employed herein, the vector V d 
(describing the velocity of m d relative to m 0 ) is a constant; and, by the simplifi- 
cations introduced, the vectors V 2 and V 3 are (essentially) directed along the 
hyperbola's asymptotes. In this regard, then, 

* (V d , V a ) + a = * (Vj, V d ), (III. 59) 

where a is defined by Eqn. (III. 57). Also, since V 3 + v d = v td (where V Xd is the 
heliocentric velocity of the space vehicle at exit from the sphere of influence) then, 


v t d = Vv| + V3 + 2V, -V, 
or 

\= V V| . V’ + 2V, v d cos (V,, V„) (III .60) 

wherein the cos (V 3 , V d ) can be obtained irom Eqn. (III. 59); and, the angle $ (V d , 
used above, has been noted to satisfy the relation 

MV d .V J ),S Hd+7Hd+ |5, d |, 


with 5qi d = - q> Hi . 

The magnitude of V Td is known, now; however, to complete a description of 
the vector its direction must be specified also. This can be done by locating the 
vector (in angular position) relative to (say) V 3 . In this regard note that (See 
Fig. m.lO). 
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Figure III. 10-Sketch of the Flight Path Inside the Sphere of Influence, with Conditions at Entry and Exit Described 



or 



sin(V V d ) = |V 3 ] sin (V,, V„) 

d 




sin 


1 V 3 I sin('V 3 , V d ) 



(III. 61a) 


where $ (V 3 , V d ) is described by Eqn. (III.59). Also, as a check, this angle is 
described by 


* (V T , V d ) = cos" 1 

D 


|V„I + |v 3 | cos (V 3 , v d ) 



(III.61b) 


According to the construction shown on the figure there is an angle equivalence 
statement which is useful here; namely. 


* ( V T_* + 8 <p • (in.61c) 

D D 

where Eqns. (HI. 61) describe the angle (V T , V d ). The angles 7t d and 8 (Px define, 
(first) the elevation angle for V Td (as it relates to the "final” flight path (T)); and, 
second, ( 8 ^) describes the absolute angular separation between the two position 
vectors r d and r D . 

In describing these angles it is necessary to obtain the vector magnitude |? D | . 
A definition of ~r D is obtained (see the figure produced here) from 


r n = 


+ Pt, 


leading directly to 


r D = [T l +^D + 2r d ' Pd ]W2 ' 

hence 

r D 3 V r d + Pi + 2r d Pd cos < f d* ^d> (m.62) 

wherein $ (r., p D ) = n /2 - p (See Fig. III. 10). However, 77 is defined from a 
description of the total angle about m d ; i.e., 
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so 


277 = 2 I CPo . | + | 0 1 + tt /2 + T?, 

'tun 


^ = T " <2 1 <P '€im 1 + |5|) • 

As a consequence of (III. 63) one notes that 


(III. 63) 


* ( f d * = J- r > = 2 ^l i J + 1*1 “ 7r - 


(m.64) 


Next, as a means of describing the angle Scp T , one can make use of the following 
geometric identity; that is, 

lp D ! COS T) = |r D | sin Sqp T , 


thus 


S<P T = sin ’ 1 


|p D l COS 17 



(in.65) 


(Note that the sign of Scp T has been arbitrarily set by 0 < 8 <p T <tt/ 2, then Scp t = 
|<p D - cp d | , where the $ <p. locate "D" and ”m d ” relative to some consistent peri- 
centric origin). 


Knowing 77 (Eqn. (m.63)) and Scp T (Eqn. (III.65)) then the elevation angle 
(y T ^) can be obtained from Eqn. (m.61c). This angle locates the velocity vector 
(V^) relative to the local horizon k 3 it is described for the final trajectory arc, 
T . Recall that this final arc extends from "D" to the terminus (which is located 
at m f , and/or its time position in heliocentric space). 

The one quantity left to be described - for the trajectory through the sphere 
of influence - is the time cl flight. This will be deferred for the moment, but 
will be discussed subsequently. 


m.E.lO. THE TERMINAL TRAJECTORY 

The terminal trajectory extends from point "D”, on the imaginary boundary 
of the sphere of influence, to the orbit of the terminal body (m f ). Because of the 
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encounter with the massive body (m d ) this final arc is not a continuation of the 
pre-encounter arc (the Hohmann ellipse). Actually the terminal path has a dif- 
ferent eccentricity and size, hence the latter orbital parameters (e T , a T ) must 
be determined. 

In regard to these parameters it should be apparent that there is (now) enough 
information known about the terminal arc to specify values for each. For instance, 
employing the expression for eccentricity 


2E. h 2 

€ = /I +— - — , 


and recalling that this flight path is described as having constant and h — in 
heliocentric space - then it can be shown that 


€ T 


r V 2 

r D V T, 




£ _ 1 I cos 2 y_ + sin 2 y_ 

' t d *d 


1/2 


(III. 66) 


Also, from the energy equation it is easy to show that the trajectory’s charac- 
teristic length is described by 


a r = 




(III. 6 7) 


2 


M 0 /r D 

m.E.ll ANGULARITY OF THE LINE OF APSIDES 

As a next task, to find the amount which the terminal apsidal axis has been 
shifted, relative to the initial (Hohmann) figure, the following method is suggested. 

By means of the conic equation, written for the terminal arc, the appropriate 
position angle (<p T ) for the radius ? D is obtained as: 


qp_ = CO s 
t d 


- 1 


el - 1 


L r D e T 


1 

e_ 


(III. 6 8) 


f 
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Now, define the angularity (change) for the apsidal axes as, Acp; which is the 
angular separation between the pericenter locations for the two trajectories ( ) h 
and ( ) T ! A sketch illustrating this shift in position(s) is shown below. 


EXTRAPOLATED PATH (T) 



Figure 111.11— Sketch Illustrating the Angular Positioning of Pericenter (P) 
on the Original (Hohmann) Path, and the Final, Post-Encounter Trajectory. 


In agreement with the figure shown above the angle locates r D relative 
to the pericenter (P T ). This, in turn, is separated from the original pericenter 
location (Pj) by the angle Acp ; and, in addition, one notes that 

op Td + A qp = cp H ^ + 8 <p d + S qp T , 

hence 


Acp = (cp H ^ + §<p d + S*P T ) - <P Td > (III. 69) 

where these various angles have been determined previously. One should note 
that Acp > 0 (dependent on the angle summation noted above). What is shown on 
the sketch (above) is a 'supposed* positive valued Acp! 
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The last parameter needed to complete the trajectory analysis is a description 
of the position on the terminal path "T" where the spacecraft will intersect the 
orbit of m f ! Knowing the orbital radius, < r f ), of this planet (or other celestial 
body), as an a priori, then by means of the conic equation it is easy to define the 
appropriate position angle as 


cp T = cos' 
T f 


V 6 ? - 1 

r f e T 


(in.70) 


where the absolute value is employed in the event that the final trajectory may be 
either hyperbolic or elliptic! 


m.E.12. TOTAL TIME OF FLIGHT 

In calculating the total time of flight it should be recognized that there are 
three (3) arcs to be considered; hence, three equations for time must be included 
in this evaluation. Since one or more of these arcs may be for either an elliptic 
or a hyperbolic arc then the two appropriate general expressions for time are: 

(a) For an elliptic arc (time measured from pericenter passage) ; 


,3 - e /| e 2 - 1 sin cp 

L 1 + € COS Cp 


+ 2 tan" 1 


1 - e . cp 

tan _ 

1 + e 2 


(HI. 71a) 


(b) For the hyperbolic arc (time measured from pericenter passage); 


* _ / a 3 €/|e--l. sin<p / e - 1 o> 

t - J : 2 tanh' 1 ( -J z tan 

> /it l + e cos cp \ f 1 + e 2 


■)]• 


(III. 71b) 


/111 — ~*L sin 9 In ( . {Z± 1 + ZEU • tan(<P/2) 

If M l+e cos cp \/T+T - / £ - 1 • tan(cp/2) 


(III.71c) 
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It is noted that, symbolically, these formulas define time as a function of the 
orbital parameters and other quantities, as 

t = t (/z, a, <p, e). (III. 72) 

The three time intervals, making up a complete description of the mission 
flight time, are obtained as the time required to fly the three separate arcs 
making up the complete mission path. Of course, this statement is hi agreement 
with the assumptions set forth earlier. These three arcs are: 

(1) The Hohmann path length, from launch at t% to intercept with the sphere 
of influence (at I); this considers an arc from a heliocentric position <p = 0 to one 
denoted as cp = cp u ! 

H i 

(2) The hyperbolic arc, within the sphere of influence, is that arc whose 
angular extent is from a planetocentric position of - <p^. to one denoted by + cp^ im ; 
and, 


(3) the arc along the terminal trajectory, extending from the exit heliocentric 
position, <p = <P T , to the terminus position, cp = % . 

D f 

In order to describe these time intervals from the expressions above the 
computations are carried cut in the following manner: 

(a) For the first arc (along the Hohmann path) the time is described via 
Eqn. (III. 71a) where the quantities employed in that expression are noted to be 

* = t (/x© V <p HI , € H ) . (III. 73a) 

(b) Within the sphere of influence the arc is a hyperbolic section, thus use 
is made of either Eqn. (III. 71b) or (m.71c). Also since this arc is symmetric 
about the pericenter (and the figure’s axis) then the time of flight would be twice 
the value calculated from either of the expressions noted above. Thus, 

♦ 

t = 2t( j u d , a, e); (III.73b) 

wherein, in reality, cp^. is the angle corresponding to a radius p l in planetocentric 
space. 

(c) For the final arc the trajectory may be either an ellipse or hyperbola - 
this can be ascertained from Eqn. (III. 66), where the eccentricity is described. 

As a consequence, this time must be calculated from the appropriate expression 
given by Eqns. (III. 71). Also , since the arc covered by the spacecraft on this leg 
of the mission is described by 
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CP- < cp < cp_ 

T D - T f 

then it will be necessary to compute the flight time by the scheme noted below: 
The time increment of interest here is, 

t = At = t f - t D , (m.73c) 


wherein 


tf = t (m©> £ T^’ 

and 

tp = t (M©* ®x* ^Tjj’ e T^‘ 

ni.E.13. ENERGY CHANGE DUE TO THE SWING-BY MANEUVER 

Calculations carried out in the sections prior to this one have been concerned, 
primarily, with the physical and geometric aspects of this overall maneuver. 
However, one of the more important parameters associated with this operation is 
still to be determined - that is, the energy change which has been brought about 
by the "swing-by" itself. 

For a determination of the energy change one should recognize, first, that 
the description is concerned with the energy change referred to heliocentric 
space; and, also, that the considerations are in regard to both the kinetic and 
potential energies referred to this space. 

Defining the change in energy, symbolically, as 

AE,=E, -E, (III. 74) 

where, typically, for example, 

E -JS-. — .i— (in. 75a) 

T 2 r D 2a r 

(again, the ± have been written to account for the possibility of either an elliptic 
or hyperbolic terminal pate); and where, for example, 
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(Ill . 75b) 


, V H d M 0 _ 

'’ H 2 "r, " 2a„ 

(written, as such, since the initial path was described as a Hohmann segment;. 

From Eqns (III. 75) one can define the change in both kinetic energy and potential 
energy (per unit of spacecraft mass); or, alternately, one can define the change 
in total energy, per se, by simply describing the change in (i = T , H ). For either 
scheme chosen a consistent definition of AE t would be obtained. 

To some degree the effectiveness of the hyperbolic encounter can be measured 
in terms of the heliocentric speed gained for a trajectory passage close to the 
massive disturbing body. In this regard, an indication of such is given by com- 
paring the quantities V H<) (pre-encounter speed) and V Td (post-encounter speed). 
Also, for that case where the final path is highly energetic (E 1t > 0) it is of some 
value to note the magnitude of the hyperbolic excess speed along the terminal tra- 
jectory. This topic will be briefly mentioned below. 

Recalling that hyperbolic excess speed is defined as 



then for the terminal path this quantity would be 



(III. 76a) 


or, as related to escape speed for this trajectory (referred to some convenient 
position) , 


V 2 = V 2 - 2 Hi. . V? - v 2 , c (m.76b) 

T t d r D T D 

since V L =2 Mo /r (at point "D", in heliocentric space). Necessarily if, from 
Eqn. (III. 76b), the resultant is negative then the vehicle did not attain hyperbolic 
speed. On the other hand, if the resultant is positive then the escape, and hyper- 
bolic speed, was reached by the vehicle during the encounter, and as a consequence 
of it. 
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IV. REMARKS 


In the foregoing sections two approximate schemes, which allow the analyst 
to study a Swing-By operation, have been outlined. Of these, the Mode I method 
i3 simplest to use and manipulate; while Mode n is a more exact procedure inso- 
far as the basic assumptions used herein are concerned. 

Needless to say these methods are not "exact" in a physical and mathematical 
sense; and, as such, both schemes will suffer some loss in accuracy under certain 
conditions. To illustrate that the methods used herein are not in full agreement 
with the physical situation; recall that the assumption has been made whereby 
m d is presumed to be fixed in a spatial position when the Swing-By occurs. In 
many instances such an approximation may not seriously affect the result; yet, 
in other cuses this assumption n ight lead to unacceptable errors. One means 
of correcting for this (one) objection would be to properly account for the 
rotation of the m d velocity vector, V d , at exit from the sphere of influence. Know- 
ing the time of motion in planetocentric space it is immediately evident what the 
magnitude of this angularity would be - and, of course, its direction is known. 

Such a correction may (or may not) be of consequence; this will be a matter for 
the analyst to decide. 

In general, the methods presented here have been developed to aid in the 
evaluation of planetary swing-by maneuvers, as a means of performing long 
distance-and-time missions. It is not suggested that these schemes will provide 
the accuracy needed for (say) guidance and precise performance calculations; 
however, the various influences which do arise as a consequence of the "Swing- By" 
can be determined, and evaluated in a not-too-rough manner. 

Some of the more subtle aspects of this maneuver (type) can be found in the 
bibliographical 1< stings; the interested reader is urged to consult these for this 
information. 
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VI. SYMBOLS 


a semi-major length of a conic (trajectory) 
b semi-major length of a conic (trajectory) 

Ej specific energy of a particle in motion about a primary 
AEj a change in specific energy 
F location of a conic's focus 
G universal gravitational constant 
h specific moment of momentum ( = | T * V | ) 
i, j, k indices 
m mass 

m d mass of the disturbing planet (about which "swing-by" occurs) 

m solar mass [ See Eqn. (111.9) ] 

© 

p orbital parameter (semi-lalus rectum) 

P period of motion; corresponding to given values of a, p. 
r radial distance (radius) 

R a radius [ See Eqn. (m.9) ] 
r IN radius of the sphere of influence 
t time 
V speed 

V 2 .(V 2 ); heliocentric velocity (speed) at entry to; and exit from; 

V,,(V_) toe sphere of influence 

3 

V H ,V T heliocentric speed, at m d , cm the Hohmann (H) and Terminal (T) tra- 
d d jectory arcs 

a relative inclination between both asymptotes of a hyperbolic path 
J3 angular position of conic's asymptote relative to limit radius (p t ) 
y elevation angle for velocity 
T speed ratio (See Eqn. A.9, Appendix A) 

A( ), S( ) increment values ( ) 

8 jd angularity of Vj (i =2,3) relative to V jd ( j= H , T) 
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e eccentricity 

6 angle between radii (p , 7 d ) 

M gravitational parameter (^ Gm; - Gm Q , = Gm s , p d = Gm d ) 

£ asymptote positioning angle (See Appendix B) 
p,(p) planetocentric radius, measured relative to m d as primary 
q> polar positioning angle (measured from pericenter) 
y angularity of asymptotes, referred to a conic's axis; see Fig. III. 3. 
co a positioning angle (Appendix A) 


Subscripts 


apo corresponding to apocenter 

D refers to exit point, leaving sphere of influence 
d refers to disturbing mass (m d ). 
f corresponds to final position (orbit, mass, etc.) 

H refers to the Hohmann arc 
h a heliocentric rc f erenced quantity 
I refers to the initial (entry) point to sphere of influence 
pe r i , P corresponding to pericenter 

oo refers to extreme (limit) position on a hyperbolic path 
( ) quantities associated with V 0 , V, 

© refers to earth quantities 
© refers to sun (solar) quantities 
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APPENDIX A 


Geometry for the Planetary Swing- By 

The planetary swing-by, per se, is considered to be confined to that region 
in space which lies within the planet’s sphere of influence. In a simple definition 
of this region, the radius of this sphere is defined by 


r 


IN 


= R 



(See Eqn. (in. 9)), wherein: 

r 1N = radius of the sphere measured from m d 
R = heliocentric radius locating m d (wrt m 0 ) 

m d , m Q = mass of the planet, and sun, respectively. 

The presumption made here is that the swing-by trajectory is a hyperbola; 
and, consequent to the definition above, the conic will have geometric and dynamic 
symmetry; with the planet, m d , (hence ti = Gm d /r 2 ) serving as the focus. 

For the geometric and kinematic conditions to be described next, concern is 
centered within the sphere of influence; though conditions at the boundary are 
also of consequence. In this regard the velocity vectors (see Fig. A.l) at the 
entry (I) to the sphere, and exit (D) from it are of interest. Also, these vectors, 
as referred to heliocentric space and planetocentric space, are of equal im- 
portance. For these purposes then, the vectors 

V. , V . (i = I, D; etc) 

h d 

are defined, wherein ( ) h refers to heliocentric space, and ( ) d implies quantities 
referred to the disturbing mass (m d ). 

Since the planet (m d ) is in motion, in heliocentric space, then its velocity V d 
is described and used here. For purposes of this analysis I V d I is constant; and 
one of the assumptions to be employed here will suggest that V d is also constant, 
but during the time of the planetary encounter only. 
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Figure A. 1— Geometry of the Hyperbolic Path Inside the Sphere of Influence 
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Certain angles, seen in Fig. A.l, are of interest here; those of primary con- 
cern are listed and described below: 

a = turning angle for the velocity vector during passage through the sphere 
of influence; 

\y = inclination of the asymptotes relative to the axi? of the planet centered 
hyperbolic trajectory; 

co = the angle of inclination for the perpendicular to the asymptote, measured 
from the figure’s axis of symmetry. (Note tha* co + \y = tt/2.) 

A satellite approaches the sphere of influence, along its heliocentric flight 
path ( approach trajectory ), and intersects this imaginary boundary at I . The 
heliocentric velocity at this point is V Ih , while the planetccentric velocity is V Id . 
These velocities are related according to 

\ = \ + V d . (A.l) 

Looking at Figure A.l it is noted that several of the angles which are of 
interest here are related by the following 


i 


\ 



and 


hence 


a + 2 4* = tt. 



a = 2 co . 


Also, it is seen that 


tan ? = b/a, 


while 


tan co - a/b, 

I 

§ 

I 


(A. 2) 


(A. 3) 

j 

} 

j 

i 

ff 

I 

i 
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hence 


t an co = cot 

which is also evident from the second relation given in Equation (A. 2). 

Since a = 2co , then from Equation (A.2) and the definition of b it follows that 


tan a. ■= tan 2oo = 


2tan oo 
1 - tan 2 co 


2 F 

1 - — 
b 2 


/? ~ 


e 2 - 2 


(A.4) 


also, with 


a 1 

tan &) = — = 


b " r ?~\ • 


(A.5) 


then 


tan V = - = /e 2 - 1 . 
a 


(A. 6) 


From the conic equation, it is noted that the radius to pericenter, defined as 
r(q> = 0) is. 


rp = — s a(e - 1); 
1 + e 


similarly, since the semi-minor length b = a/e 2 - 1, and the eccentricity can 
be determined from e 2 = 1 + b 2 /a 2 , then it follows that 





or, after manipulation 


b ^ 



a . 


r P 


(A. 7) 
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Equation (A. 7) is important to these problems since it is (both) the semi-minor 
axis length as well as the perpendicular distance to the asymptotes from the locus, 
F, (see Fig. A.l). 

The semi-major axis length (a) is also significant in that according to previous 
assumptions the vectors V ld and V Dd are essentially parallel to the asymptotes. 
Consequently these play the roles (approximately) of the so-called "hyperbolic 
excess velocity" vectors (V®), which are defined as being parallel to the asymptotes 
at the extremities of the hyperbola (as r - oo). 

One should recall that this velocity vector has a magnitude defined by 



(A. 8) 


Next, a parameter (r , describing a speed ratio, will be introduced: this 
quantity will be defined (h me) as 


r* § _XL , 

M/tp 


(A.9) 


where /u / rp is the circular speed for an orbit of radius rp (the pericentric 
radius). As a consequence of this definition it can be shown that the following 
expressions occur: 


(1) Using Eqns. (A.8, A.9), 

V ? 

-L- - r- 2 ; where H 2 = — L. (A.IOa)* 

r p d d /V rp 


(2) From Eqns. (A.7, A.9), 


b 



or, b = 



(A. 10b) 


* 


One should note that 




, since 
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(3) From Eqns. (A.5, A. 9), 


tan ¥ = cot w = Pj ^2 + Tj 2 ' . 

d d 


(A. 10c) 


(4) Combining Eqns. (A.8 t A. 9, and A.4), 


tan a = 


r .p 2 + r .p- 1 


and 


(5) Utilizing the Eqns. for {rp) and (A. 8), 

« = 1 + r* 


(A.IOd) 


(A.IOe) 


Note: these expressions apply only for the hyperbolic encounter. 

From these expressions it is apparent that for a fixed value of rp , as V Id is 
increased the angle a decreases. That is, the "turning" of the velocity vector 
is reduced; this can be attributed to a reduction in time during which the space- 
craft can be influenced by the attracting (or fly-by) mass. Therefore, to obtain 
a maximum of turning (a) it is essential that a slowest approach velocity (V^ 
should be available, at a fixed value of rp . Conversely, for | Vi d | fixed the value 
of a increases as rp is reduced; hence, a close passage to the attracting planet 
is desired if a is to be maximized (for a given value of | V, | ). 

In substance, then, the maximum turning (a) occurs when rp and IV are 
smallest, compatible with the assumptions and constraints of the mathematical 
model assumed here. 

It should be recognized that there are some least values for quantities con- 
nected with this study. For instance, the smallest value of rp which can be used 
here is that corresponding to the dimensions of the attracting planet (and its 
atmosphere, in those cases where atmospheres are significant). Secondly, the 
entire problem has pre-supposed that the trajectory, during swing-by, is hyper- 
bolic; therefore |V | - and lv Dd l - must be at least as large as the parabolic 
speed relative to the attracting planet. This assures the geometric form of the 
conic assumed earlier. 
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One final remark before leaving this section; this is in regard to the distance 
FP(= B) shown on Fig. A.l. It has been stated (without proof) that this length is 
equal to the semi-minor axis (b) of the geometric figure; actually this is easily 
demonstrated, and will not be developed here. 

Since the hyperbolic trajectory is assumed to be generated as a two-body 
flight path, it should be evident that the usual two-body relations are applicable 
and can be used to describe this geometry. Also, the state of motion - as well 
as time of flight - relative to the planetary mass (/z), can be described from 
two-body results. 


Velocity and Energy Change from the ’’Swing- By” Maneuver 


Figure (A.l) shows the heliocentric exit and entrance velocity vectors as 
and V Ih , respectively. Consequently the change in velocity (8 V), due to an 
encounter with an attracting mass can be defined as 

§V = V n -V T . (A. 11) 

u h x h 

Making use of (A.l) it is evident that 

sv =wvv <A - i2) 

as seen from the sketch here, (Fig. A.2). 



SPHERE OF 
INFLUENCE 


Figure A.2_^Sketch illustrating the velocity 
change, 8V = Vpj - V|j. Here the Circle 
of Influence is a cross-section of the Sphere 
of Influence. 
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The magnitude of this change ( 1 8V|) is obtained as 

AV = | S V | = [(V- -V, ) • (V ! 

H d do 

or, expanding and recognizing that V D = Vj , then 

d d 

AV = V, /2(T - cos a) * V D , A <T - cos a). (A. 13a) 

*d d 

An equivalent expression, more compactly written,* is 


AV = 2V. sin ± . 

! d 2 


(A. 13b)* 


The total energy, relative to the hell': center, at entry to the sphere of in- 
fluence (at I) and exit from it (at D), can be written, generally, as 

V. 2 

E. = _ if® = - it®, (i I, D'i 

* 2 r. 2a. 

l h * 


where /*© is the gravitational constant for the sun; r ih is the radius to the points 
(i) from the sun; and, "a. ” if» the te mi- major length describing the heliocentric 
trajectories (before and/or after encounter). In general the energy is not the 
same before and after the swing-by - this is normally one of the reasons for 
attempting this maneuver - therefore, a change in total specific energy, 5E i >< 0, 
is expected. 

Under the assumptions made for this analysis, it will be presumed that in 
heliocentric space, r Ih -v r^; therefore, the change in energy may be expressed 
(approxi mate ly) by 





(since the r 4 are assumed equal). 


(A. 14a) 


*lt wos noted earlier that a/2 = a>, and th at ~ /u/ a, hen ce on approximate evaluation of AV 
is obtained as AV = 2VJJ sin a/2 2 //t o sin co = 2 /JI7a cos ty. 
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In this expression. 


v 2 - I V + V i 2 - V 2 + V 2 + 2 V -V 

\ 1 »d + jl ■ ” d d D d j 


and 


V ? = |V. +V„P = vf +V2 + 2V V d ; 




so that, accounting for the angles ( j = I , D) on Figure A.2; also recalling 


that V D =? Vj then, 

d d 


also 


VS = V? + 2 v d V cosf D + » V| + 2 V V d cos + V*; 

h i a d d 


V? = V l * 2 V V rf cos + V d * % + 2V V d cos + V^i 

h d d da 


aiKl, therefore, 


SEj y v d Vj (cos £ D - cos ^j) 

d 


(A. 14b) 


which is equivalent to 


5E i - v d ' (V - V) » V d • 8 V 

a d 


(A. 14c) 


(See Eqn. (A. 12)). 

It should be evident from these expressions that SE 1 >< 0 depending on the 
orientations of V Ih and . Another means of stating this is to note that the 
sign of SE 1 is compatible with the direction cosine of SV relative to V d (See 
Eqn. '* .14c)); or, it depends on the sign and/or magnitude of cos and cos 
(See £qn. (A.14b)). In any regard, it is seen that the energy change, relative to 
the heliocenter, can be described as a gain in energy ( SE. > 0) or a loss of energy 
(8E t < 0), for the satellite, depending on the approach to and retreat from the 
sphere of influence. If SV is generally in the direction of V d then energy is 
gained ; however, if the direction of SV is in opposition to V d , then the vehicle’s 
passage will represent a giving up of energy to the attracting body, or a ’’loss" 
of energy relative to the heliocenter. 

Before leaving this section, it is desirable to express some of these results 
in terms of the defined speed ratio ( T ). In this regard, it should be recalled 
that d 
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1 


, a 

tan co - — ~ 

b 

then according to the sketch (below) 



sinw; ; 

and, according to the footnote with Eqn. (A. 13b). 



Extending this to encompass the statement given as Eqn. (A. 14c) it is evident 

that 


SE i - v d 


oV 
S V | 


2 






1 + ^ < 2 + ^ 


(A. 16) 


where (now) the magnitude of the change in velocity, and the change in energy, due 
to the encounter with an attracting mass, have been expressed in terms of known 
(or selected) quantities, and the speed ratio (C ). 
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APPENDIX B 


Extremals for the Change in Velocity and Specific Total Energy 

Having evaluated the changes produced by a "swing-by", the next logical 
step is to investigate these to ascertain the conditions and magnitudes for the 
extremals of each. Taking these in logical order, the change in velocity (magnitude) 
will be determined first. 

In order to describe an extremium for the change in velocity the magnitude 
of SV will be studied first since the direction of SV depends on the vector 
difference V D d “Vi d , wherein each of these vectors has the same magnitude. 
Generally the extremal will be dependent on Vi d ( = V Dd ) since the pericenter 
distance (r p ) is dependent on the energy (hence the entry speed) for the hyperbolic 
path. In this regard, then, the extremal for AV will be determined in terms of 
V Id (or, equivalently, in terms of r x — and, in particular, in terms of since 
this is the order of this quantity as it appears in Equation (A. 15)). 

Therefore, differentiating Eqn. (A. 15) with respect to Pj* and setting the 
result to zero leads to the expression 



M/rp 

[i + i7 d (2 + r, 2 d )P 


(l-r, 2 ) = 0; 

d 


from which the condition for the extremal is 


r x a = l, (B.i) 

l d 

and/or 

V 2 _ M . 

I . “ 7— 

d r p 

Necessarily this suggests a condition for a maximum since the minimum for 
AV is obviously zero. 

Equation (B.I) implies that if AV is to become (AV) max , then the speed(s) Vj 
(and V Dd ) should be equivalent to the circular speed corresponding to an orbit d 
of radius rp ! 
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Using Equation (B.l) in Equation (A. 15) one finds the maximum of SVto be 


AV 


max 



(B.2) 


hence the magnitude of the maximum gain in velocity is equivalent to this (circular) 
speed also. The direction of §V has been determined as the vector sum of V D d 
and -Vi d ; these (in turn) depend on the entry and exit positions at the sphere of 
influence. Note that the gain (which is depicted here) is independent of the planet's 
velocity (V d ) and is solely determined by the closeness of the passage to the planet 
(rp), as well as the size of the attracting mass (m). 

The extremal for the change in specific energy is more involved since (see Eqns. 
(A.14))there is an implied conditionality for the orientation of vectors involved here. 
Certainly the extremal for this quantity depends on both the direction and magni- 
tude of SV. For the moment, neglecting the magnitude of SV , it is evident that 
SEj has its largest value when the angle between v d and SV is zero; and, in 
contrast, SEj is a least value when these vectors are directed opposite to one 
another. Thus, the extremals of £E 1 are dependent on the statement 

cos(V d ,SV) = ± 1. 

Now, having defined a magnitude, (AV) max , then it would appear, from Eqn. 

(A. 14b), that the extremals for SEj are defined according to the conditions 

(1) AV =? (AV) max (B.3) 


and 


(2) cos £ d - cos = ± 1. 

Condition (1) has been described, see Eqn. (B.2); condition (2) can be established, 
however, it should be noted that in agreement with Figure A.l, 

= a » 

therefore, condition (2) is better stated as 

cos £ d - cos(£ d + a) = ± lj (B.4) 


67 


ms- 


or, expanding, 

cos £ d (1 - cos a) + sin £ D sin a = ± 1. 

After a value for a is decided upon, then it becomes a relatively simple matter 
to evaluate Eqn. (B.4) for the angle ; and, to use Eqn. (B.3) to ascertain a 
proper value for . 

It should be apparent that the angle a is not arbitrary; as a matter of fact 
Eqns. (A. 13b), (B.l) through (B.3) lead immediately to the result 

m m . x = (2V sini) 

\ * /max 

hence, 

HE. 2 [E sinl ; 

V r P V r P 2 

or, the requirement that is implied here is 


(B.5) 


Now, solving Eqn. (B.4), accounting for the proper value of a , it is found that 


and, therefore, 


cos VJ 

2 + 2 


sin £ d = ± l 



Now, using Equation (B.3), it is observed that the corresponding values of 
£ , needed to extremize SEj , are obtained from 

1 

cos = T j 


and, as a summary, the following table is presented. 


Table BJ 

Conditions for Extremizing SEj Due to a Planetary Swing-By Maneuver 


SEi 

cos (v d , 8V) 


Comment 


0 

7t/3 

277/3 

Energy is added for the 





satellite's heliocentric path 


7 T 

277/3 

77/3 

Energy is extracted from the 





satellite's heliocentric path 


Graphically, the solutions noted in the above tabulation correspond to the 
following fly-by conditions, in the sphere of influence: 




(energy is. added to the (energy is extracted from the 

heliocentric path) heliocentric path) 

The information above indicates that the extremals for SEj are obtained by: 
(1) a gain in energy by passing behind the attracting mass; and doing so by keep- 
ing the vectors V d and SV parallel (directly), also having AV = (AV) mtjx . The 



second extremal, (8E 1 ) min , occurs when the satellite vehicle (or attracted mass) 
passes in front of the attracting mass; this occurs when 8V and V d are oppositely 
parallel, but while AV is retained at the level of (AV) max . All of this is, of course, 
based on the assumption that ip has been pre-selected, a priori, and that this 
parameter does not otherwise influence the extremals. 

Apparently conditions relative to the direction of Vlh ■ V , etc., the size of 
the trajectory flown during th e pas sage through the sphere of influence, and the 
magnitude of V d (other than vV Ap ) will produce changes in 8V and 8Ej but not 
the values corresponding to the extremals unless the necessary conditions for the 
extremals are satisfied. 


